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COSMOLOGY FROM THE 
GENERALIZED EINSTEIN EQUATIONS 

Recently, we have investigated cosmological applica- 
tions and confronted them against astrophysical data the 
following class of gravitational Lagrangians: 



of ACDM model is described by V\cdm = — \ (Aa 2 + 
rjcf '' 
ter. 



) where tj is a density parameter for the dust mat- 



As it was shown in [1] the equation (2) leads to two 
classes of cosmological models implemented by differ- 
ent solutions of generalized Einstein equations. 



L = Jg(f(R)+F{R)L d ) +L mat = 
s fg(R + aR 1 + PR [+S +R [+a L d ) +L ma , 
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within the first-order Palatini formalism [ 1 ] . Here = 
— ifg^ dn^dyty is a scalar (dilaton-like) field Lagrangian 
non-minimally coupled to the curvature and L mat repre- 
sents perfect fluid Lagrangian for a dust (non-relativistic) 
matter. The numerical parameters a , j3 , 8 , (7 are to be de- 
termined by astrophysical data. 

Applying (Palatini) variational principle compiled 
with flat FLRW metric one arrives to general Friedmann 
equation: 



2(f + F'L d ) [3f - f'R + (3F - F'R)L d ] 
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where H = - denotes the Hubble parameter related to the 
FLWR cosmic scale factor. This reconstructs the ACDM 
model under the choice f = R — 2A, F = 0, which is the 
limit a = 0, 5 = - 1, j3 = 2A. Setting further A = leads 
to Einstein-de Sitter (decelerating) universe. 

We want to recall that the generalized Friedmann 
equation under the form: 
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(which is always the case for the Palatini formalism) 
leads to one-dimensional particle like Newton-type dy- 
namics which is fully described by the effective poten- 
tial V(a) — — ^a 2 G(a). This relevant property allows us 
to compare various cosmological models on the level 
at the effective potential functions and the correspond- 
ing phase-space diagrams. Particularly, the dynamics 



Model I 

Solving equations of motion by 

R = p = na~ 3 , a = -8 



(4) 



one obtains generalized Friedmann equation under the 
form 
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where 



(6) 



are dimensionless (density like) parameters. 

Model II 

Another cosmological model can be determined by 
'7 
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.(1-8)1) 
which leads to 
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where now 
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£lo,a = aHgQop (9) 



Both models have CIq m ,Q.Q p,8 as free param- 
eters. By the normalization condition H(0) = Hq, only 
three of them are independent (Hq denotes the Hubble 
constant). 



FITTING PARAMETERS OF THE 
MODELS 

In order to estimate the parameters of our models we 
use a sample of N = 557 supemovae (SNIa) data [3], 
the observational H(z) data [4], the measurements of the 
baryon acoustic oscillations (BAO) from the SDSS lumi- 
nous red galaxies [5], information from CMB [6] and, as 
an adition to [1], information coming from observations 
of GRB [2]. 

The entire likelihood function Ljot is characterized 
by: 

Ltot = LsnLh-LbaoLcmbLgrb- (10) 

We have assumed flat prior probabilities for 
all model's parameters. We also assumed that 
H Q = 74.2 [kms^Mpc- 1 ] [8]. 

The likelihood function is defined in the following 

way: 

_ ^fheor_ /i obs\2l 
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where: a, is the total measurement error, pL } 
M is the measured value (m,— apparent magnitude, M- 
absolute magnitude of SNIa), J u ! - W = 5 \og w D u + JZ = 
5 log 10 d Li + 25,Jt = -5 log 10 // + 25 and D Li = H Q d Li , 
where du is the luminosity distance given by du = (1 + 
Zi)c Jq (with the assumption k = 0). In this paper 
the likelihood as a function independent of Hq has been 
used (which is obtained after analytical marginalization 
of formula (11) over Hq). 

For the H{z) data the likelihood function is given by: 



= exp 
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where H (zi) is the Hubble function, //, denotes observa- 
tional data. 

For BAO A parameter data the likelihood function is 
characterized by: 



l B ao ex P 
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where A' 1 
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0.469 ±0.017 for z A =0.35. 

We also use constraints coming from CMB tem- 
perature power spectrum, ie. CMB R shift parameter 



[7], which is related to the angular diameter distance 
(Da(z*)) to the last scattering surface: 



\/Cl m Ho . . . 

R= (1 +z*)D A (z t ). 



The likelihood function has the following form: 



Lcmb <* exp 



1 (R-Robs) 2 



(14) 



(15) 



where R obs = 1 .725 and a A 1 = 6825.27 for z* = 1091 .3 
[6]. 

The likelihood function for GRB data is defined as: 



Lgrb °= exp 



^- n th ( Zi: n. m ,a A , 
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The mode of joined posterior pdf as well as mean (to- 
gether with 68% credible interval) of marginalized pos- 
terior pdf were calculated, by means of Markov Chains 
Monte Carlo analysis, using free accessible CosmoNest 
code [9] which has been modified for our purpose. The 
results are presented on fig. 2,3. 

The numerical values of best fitted parameters for 
two our models as well as for ACDM are collected in 
table 1 : the previous estimations without the GRB data 
(i.e. SNIa, H(z) and BAO and CMB) are shown in top 
part of the table. The new estimations including the GRB 
data occupy bottom part of the table. 

Quality of the estimation can be visualized on the 
Hubble's diagram (fig. 1). Both of our models are in good 
agreement in the observational data. 




FIGURE 1. Comparison of Hubble's diagrams for models: I 
(blue), II (magenta) and ACDM (black). 



CONCLUSIONS 

In this paper we continued and completed analysis of 
new cosmological models which were previously de- 
scribed and investigated in our paper [1]. Adding GRB 
data [2] allowed us to obtain better constraints of param- 
eter £2 a which wasn't present previously. 

As it can be seen on the potential plots (fig. 4,5, both 
models dynamically mimics ACDM model from the Big 



-50 50 




FIGURE 2. Constraints of the parameters of model /. In 2D 
plots solid lines are the 68% and 95% confidence intervals 
from the marginalized probabilities. The colors describe the 
mean likelihood of the sample. In ID plots solid lines denote 
marginalized probabilities of the sample, dotted lines are mean 
likelihood. For numerical results see Table 1. 
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FIGURE 4. The diagram of the effective potential in 
particle-like representation of cosmic dynamics for model I 
versus ACDM model. Note that till the present epoch two 
potential plots almost coincide. Particulary, one can observe 
decelerating BB era. Maximum of the potential function cor- 
responds to Einstein's unstable static solution. Discrepancies 
become important in the future time: e.g. discontinuities of the 
potential functions (vertical, red line) denote that V — > — °°, i.e. 
d — > <*> for a — > af'"" 1 . It turns out to be finite-time (sudden) 
singularity. In any case the shadowed region below the graph is 
forbidden for the motion. 




FIGURE 5. The diagram of the effective potential in particle 
like representation of cosmic dynamics for the model II a =o 
versus ACDM model. Maximum of the potential function cor- 
responds to unstable static solution. Again, until the present 
epoch there is no striking differences between plots. One can 
observe finite-size sudden singularity in the near future (ver- 
tical, red line). In any case the shadowed region below the 
potential is forbidden for the motion. 



FIGURE 3. Constraints of the parameters of model //. The 
meaning of the colors and the lines this same as in the picture 
2. For numerical results see Table 1. 



Bang singularity until the present time. Discrepancies 
will appear in the near future. Both of our models predict 
the final finite size and finite time singularities (at a = 
1.673 for the model I, and at a = 1.559 for the model 
II). However, comparing with our previous simulations, 
adding new GRB data has changed properties of the 
model II (Big Bounce is now replaced by Big Bang). 
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